A non-equilibrium system in a steady state: wind waves in the open ocean 
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We derive scaling laws for the steady spectrum of wind excited waves, assuming two inviscid fluids 
(air and water) and no surface tension, an approximation valid at large speeds. In this limit there 
exists an unique (small) dimensionless parameter e , the ratio of the mass densities of the two fluids, 
air and water , independently of the wind speed. The smallness of e allows to derive some important 
average properties of the wave system. The average square slope of the waves is of order |Zn(e 2 )| , 
a small but not very small quantity. This supports the often used assumption of small nonlinearity 
in the wave-wave interaction. We introduce an equation to be satisfied by the two-point correlation 
of the height fluctuations. 



I. INTRODUCTION 



Wave turbulence makes a problem interesting both from the point of view of fundamental science and because it has 
obvious connections with real life phenomena which play a significant role in many human endeavors. Wind-excited 
— i" waves provide an example of a non-equilibrium system which may be in a steady state. After many scientists, we 
consider below the following problem: a constant wind blowing on the horizontal surface of an infinite Ocean (in 
the three dimensions), and exciting, after a transient, a system of random fluctuating waves that reach a turbulent 
steady state. After Pierson and Moskowitz [l| this makes what is called a "fully developed sea" in the oceanographic 
literature. One studies the statistical properties like the variance of the surface slope as well as the strength of the 
impact of the waves on various structures and finally the power spectrum of the surface elevation. This has been 
^ . thought about for many years 0,0, mainly by using the weak interaction approximation, based on the (unexplained 
to the best of our knowledge) observation that the slope of waves is small on average. Below we approach this problem 
by using simple scaling ideas based on the existence of a small parameter, the ratio of the mass density of the air 
Pair to the mass density of water p w . Let e = be this small ratio, about 10~ 3 . Assuming incompressibility, (and 
'■^ . neglecting, for the moment, surface tension and viscosity), the data with a physical dimension are the acceleration of 
gravity g and the (uniform) wind speed U. As shown by Newton, one can make out of these two quantities a length, 
^ — g an d a time, r = jj = — . Here we shall use A and r as units of space and time, that leads to a system of 
equations without the physical parameters g, U This is in contrast with "ordinary" turbulence where length scales 
are fixed in two ways: the large scale is determined by the geometry of the flow, the diameter of a pipe for instance, 
although the small scale is the Kolmogorov scale where viscosity becomes important. Moreover we shall argue that, 
in wave turbulence, there is no step by step cascade from large to small scale where dissipation takes place. 

Consider the fluctuations of the surface elevation, 5h(x, t), which depends on space, i.e. on the horizontal coordinates 
, x (boldface are for vectors in the Euclidean geometrical space) and on time t. From the data, the only scaling 
parameters are the length A and the time r. Given e, the fluctuations 5h(x.,t) may then be written as 

o: 



6h(x,t) = \6H£,i). (1) 

> . A T 

Therefore we expect that Sh(x, t) is of order A, and changes with respect to space and time with the typical scales A, 
and t, eventually times factors depending on e only. More precisely we expect that the horizontal dependence is scaled 
with A without factor depending on e, at least for e small, because the typical wave-length of the unstable fluctuations 
is of order A. This is in qualitative agreement with the observation that, the larger the wind speed the larger is the 
average wave-length of the wave system, as reported for instance in [2|. Note that, if this average wave-length were 
much less than A (i.e. if the surface was flat at the scale A) , this wave system would be unstable against fluctuations 
at wavelength of order A, which is impossible if the wave system has reached a statistically steady state. Moreover 
this kind of assumption is fully consistent with the Pierson-Moscowitz spectrum where the wave period is of order r 
although the amplitude is multiplied by a small numerical prefactor of order 10 -2 explained below by a non-trivial 
dependence of the average wave-amplitude on e [5]- Indeed the wave-height should scale like /(e)A, where the /(.) 
function tends to zero as its argument tends to zero ( if /(0) would not be zero, the water-wave would keep a finite 
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amplitude in absence of wind, which is against common sense, because there are no waves if no instability feeds 
the wave system). Therefore the "scaled" stochastic function SH(X,T), with capital letters for scaled quantities, is 
proportional to the reducing factor /(e). We shall prove that /(e) is logarithmic with respect to e. Assuming invariance 
under translation in space and time , the pair correlation Ssh(X, T;Y,T + T') =< 6H(X, T)5H(X + Y,T + T') > is 
a function of (Y, X") and of e only. This does not mean that it is isotropic, as one expects the direction of the wind 
to induce such anistropy. 

Let us now comment on the neglect of the viscosity of air and of water, and of the surface tension a. The latter 



introduces a new length scale, the capillary length A c = i/- 2 — • Although surface tension of real sea- water is a highly 



variable quantity, the capillary length is somewhere between one centimeter and one millimeter, far smaller than 
the length scales we are concerned with. This does not mean however that its effects are negligible: surface tension 
could regularize the instabilities at the shortest scales, which otherwise could lead to the development of cone-like 
singularities on the free surface @, while white caps are observed instead. Here we are interested on large scale 
phenomena only, then capillarity effects are not considered, and the capillary length scale will be taken as zero. Here 
we consider high wind velocity regime, more precisely the regime above the onset of wave-breaking [2| which occurs 
for Uio around 7m/ s, ( Uio is the wind speed measured 10 meters above the sea surface), where gravity becomes 
the only governing parameter of wind-wave interaction. At inferior velocities, surface tension may either dominate 
gravity everywhere or just make impossible the steepening leading to wave breaking [8j. Similar statements could be 
made abouth the effect of viscosity: if the Reynolds number is large, viscosity is relevant at very small scales that can 
be taken as simply zero. 

Now we would like to discuss the possible relationship between the ideas introduced above and the one of cascade 
of energy. The idea of cascade of energy is prevalent in fully developed turbulence, where it is assumed that energy 
goes step by step from large to small scales, energy being injected at the largest scales and dissipated at the tiny 
Kolmogorov length set by viscosity. The application of this idea to wind-driven waves meets several difficulties. The 
first one is that no scale is fixed from outside by the geometry of the flow, therefore it is hard to claim that energy 
is "injected" at large scales only. The wind acts on each scale, as illustrated in the different models using an input 
source term that is linear or quasi-linear with respect to the wave-spectrum. One must note that the Kelvin- Helmholtz 
instability is effective from A, defined above, down to the capillary length A c , the smaller length scales being the most 
unstable. In such a fully nonlinear problem, the process of direct cascade (from large to small scales) has no serious 
support. Furthermore the numerical solution of the equation of weak turbulence shows instead a drift of the peak 
spectrum towards larger and larger scales 

The reverse statement, that is that wave energy is continuously injected at small scales and cascade up to large 
scales is even more meaningless: at large scales viscosity is negligible and there is no dissipation (at least if "large 
scale" is meant in its general acception) and the energy has no way to go to in a steady statistical state (this does 
not apply to an initial value problem without sustained energy input). Moreover an inverse cascade is in complete 
contradiction with casual observation: wind waves tend to propagate along the wind direction, although, resulting 
from an inverse cascade, their spectrum should become isotropic after a few steps in this cascade. To conclude this 
discussion, a step by step cascade in any direction is not the dominant process for energy transfer across different 
length scales and cannot explain dissipation. 

Another remark may be also relevant: the air flow above the sea surface is highly turbulent and one expects the 
formation of a Prandtl logarithmic layer there. The scaling parameter of such a Prandtl layer is a flux of horizontal 
momentum per unit area of the surface. Compared to our scaling via the wind speed, this introduces logarithmic 
corrections that are, most likely, hardly detectable. However this scaling via the flux of horizontal momentum, as 
compared to scaling via the wind speed, is not completely without consequences. Actually this flux should be the 
same in the water below the surface, with an horizontal momentum of the same order of magnitude as the aerial 
flux. The flux of horizontal momentum scales like pu 2 , therefore the average horizontal speed underwater is of order 
of Uy/e (U wind speed), for example a wind speed of 60 km/h should generate an underwater current of about 1.9 
km/h. This prediction fully agrees with the observations, u sea /Uia = 0.032 [10], and the common "3 per cent rule", 
i.e. that the surface current is approximately 3 per cent of the wind speed. 

Having excluded dissipation by a multi-step cascade of energy as the dominant mechanism of dissipation, we are 
led to look at another mechanism, wave-breaking. This is known to be a strongly dissipative and fully nonlinear 
process, Therefore it could look impossible a priori to estimate the power lost per unit area by wave-breaking, given 
for instance the pair correlation of the fluctuations of height. This pair correlation assumes implicitely a smooth 
wave propagation, namely a singly defined h(x 1 y,t) and a non self-crossing surface, although wave-breaking is a 
complex nonlinear process, requiring to account for physical effects like surface tension and viscosity. Nevertheless 
wave-breaking of dominant waves (whose wavelengths are near the peak of the spectrum) should be a relatively rare 
event, as we shall explain. Let us precise that we may also include in the breaking process the micro-breakers, or small 
white caps, whose density increases with the the wind speed, and saturates at Uio around 20m/s. In this regime, 
each wave of any wavelength has a white cap on its crest. Nevertheless, in the stationary state, the crest length of 
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breaking waves is a small fraction of the total crest length. Those white caps may replace the formation of conical or 
wedge- like singularity [a] of the inviscid equations. Note that white caps appear also after the breaking of a very long 
wave. Here we focus on the very large energy dissipation process due to wave-breaking of all waves. 

Let us estimate the probability of such a process. In a single wave-breaking event, the lost energy Whr scales as 

U 2 X 2 ash, i-e. the energy corresponding to a "typical" volume, A ash of a wave having an horizontal surface A , a 
height a$h, and propagating at velocity U (which is around the velocity reached by long waves, the peak waves, before 
breaking). By comparison, the kinetic energy of the wind which occupies the same volume, is smaller by a factor e. 
The balance of energy may be obtained in the steady regime when the frequency of occurence of weave-breaking is of 
order e, since in this case the power lost per unit area is of the same order, eW4,-/A 2 T. 

The result of these considerations is that, at a given time, the area covered by wave-breaking events is of order e 
times the total area, independently on the wind speed. Note that this is not the proportion of area covered by foam, 
much bigger than the area covered by wave-breaking events in high wind conditions, because of the spreading of the 
foam by the wind, and because of the white caps. 

Let us now estimate the e-dependence of the magnitude of wave fluctuations. Given the statistics of the wave 
fluctuations, and because wave-breaking is a rare event, it depends on the probability of large fluctuations. Those 
large fluctuations are rare because of the smallness of e. Said otherwise, the scaling laws for the magnitude of 
the fluctuations have to be amended in order to take into account this smallness: otherwise, if the magnitude of 
the fluctuations was such that the typical wave length and the typical wave height were both A, the probability of 
occurence of wave breaking would be of order one too, although we argued it is of order e. Therefore the amplitude 
of the wave should be small compared to A, to make exceptional the nonlinear evolution towards wave-breaking. 
This implies that, predominantly, the wave system is described by the linear approximation of the wave equations. 
Therefore, at a given location, and according to an idea that seems to have been stated first by Planck for the classical 
(Rayleigh) part of the black-body spectrum, the fluctuations of the free surface are predominantly Gaussian because 
they are made of a linear superposition of waves with a continuous distribution of frequencies. The most interesting 
quantity from the point of view of wave-breaking is not so much the amplitude of the wave, but its slope. Being given 
by a linear transformation of the amplitude it has also a Gaussian distribution. Therefore the derivative ^ = h lX , 
namely the gradient of the height along x, the wind direction, has the probability distribution 

P(M = (0 U/2 6 ^ ^ 

where a 2 , is the variance of the slope along the x direction. 

Wave-breaking is a nonlinear phenomenon where the free surface becomes first vertical and then overturns. Before 
this happens, the slope has to reach values of order one. Indeed such a phenomenon is not described by the linear 
approximation for wave propagation, since it assumes the height to be a single valued smooth function of the horizontal 
coordinates. Nevertheless we may assume that, before the wave locally overturns, its slope gets finite (non-small) 
values that are at the border of applicability of the linear approximation [|| . The probability of overturning may be 
approximated as the tail area of the density V(h iX ), assuming that the waves propagate predominantly in the wind 
direction) , 



V br = 2 / V(h x )d(h x ) = 2erfc(-=— ), (3) 

Ja 

where a is a parameter of order unity defining the limit slope (above which the wave breaking occurs with a finite 
probability), and erfc(z) = J 2 °° d<;. In the limit of large z, erfc(z) ~ ez ^^lf ^ . According to the arguments 
presented before the probability Vbr must be of order e to ensure the balance between energy input and dissipation 
by wave-breaking. In the limit e << a x << a, this gives 

(4) 



This is of order 0.27 for the case of wind-waves, if one takes a = 1. 

This behavior is quite weak, in agreement with the observation that the slope is small on average, but not very 
small [11]. In this respect the density of breaking events is far more sensitive to the smallness of e than the wave 
amplitude itself. Indeed it seems difficult to change e, so the prediction of an e dependence of the probability of 
wave-breaking is hard to test. However one may think to a value of e close to 1, with water and oil for instance (E.J. 
Wesfreid private communication), which would then yield a wave system where wave breaking is almost everywhere, 
a clearcut prediction of the present theory. 
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Let us sketch now a more quantitative approach. As has been shown since a rather long time Q , small non linearities 
yield a kinetic equation for wave turbulence. In this theory the interaction between waves of various wavenumber and 
frequencies yields a kind of Boltzmann-like theory. In the Hasselmann-expansion, written with respect to the elevation 
amplitude as small parameter, the nonlinear interaction between the waves appear as a serie involving successively 
four, six, etc.. wave- interactions 12]. Actually the relation (Q} proves that the ratio between the six and four wave- 
interaction terms is always small. Therefore using our scalings, and the small slope variance as small parameter, we 
infer that formally the Hasselmann-expansion can be continued order by order (even for non small waves), at the 
price of fastly growing complexity. There is no reason that this expansion becomes ill-defined at any finite order. 
This does not mean however that, even by including the kinetic terms at all orders, all the physics is captured. This 
is the well-known phenomenon of expansion beyond all orders Indeed the occurence of wave breaking depends 
on effects transcendentally small with respect to the expansion parameter, here (lne 2 ) -1 / 2 . The "standard" way of 
getting such transcendentally small terms is by looking at the general large order term in the expansion, getting its 
leading order part and then summing the largest terms of each order, something perhaps doable, but surely very 
cumbersome in the present case. 

Let us now include in a qualitative sense the dissipation due to wave-breaking, in the light of the above estimations 
summarized in equations ©-Q- Using the scales variables, the steady state spectrum N(K) of the surface elevation, 
which is the spatial Fourier transform of the single time correlation function Ssh(X-, T; X + Y, T)), writes 

M(k) = — f dYe iYK 5 (5ff (Y). (5) 



(2^) 2 j 

It obeys the stationary Hasselmann equation , 

S nl [Af(K)] + S m [N(K)} + S diss [Af(K)] = 0, (6) 

where the term S n i represents the exchange by nonlinear interaction between waves, Si n the input by the wind, and 
Sdiss for the dissipation by whitecaps and wave-breaking of dominant waves. 

The term S n i describing the transfer across the spectrum of fluctuations due to four wave interactions writes 

S nl [N(K)] = J dKidK^KglTo^gl^CK + K x K 2 K 3 )5(Cl + n x - fi 2 - fi 3 ) x 



A^(K 2 MK 3 MK) ^ + ^ - j^- - ^) , (7) 

where S(.) is the Dirac function. The transition matrix T0123 = T(K, Ki, K2, K3) has a rather complex explicit 
form [l4j]; and is equal to K 5 / 2 times a numerical function of the ratios K/Ki, with i = 1,2,3, and of the angles 
between the four vectors (K, Ki, K2, K3). Moreover il = Ty/gk , and Oj = T^fgki . The input from the wind is just 
eJ\f(K)Kx, where Kx — k x is the Cartesian component of (k) along the wind. This simple approximation could be 
refined by multiplication by a function G(K) representing in dimensionless notations the detailed dependence of the 
input as a function of the wave number of the waves. We just take G = 1 to make the exposition simpler. 

An essential point of our formulation is that, in our dimensionless variables this term is small, proportional to 
e, this smallness being a consequence of the scaling laws, not assumed from the beginning. Moreover this term is 
proportional to k x to represent the angular dependence of the rate of growth of the Kelvin-Hclmholtz instability, and 
it is also proportional to the intensity of the fluctuations because this is a linear instability. 

The energy loss by wave-breaking is equal to 6'7 7 6 r A/'(K) because it is proportional to the small probability of this 
process, and it should be proportional to the spectrum itself. The numerical factor b' is of order 1, as discussed below. 
Using the expression ([3]), with er 2 = f dKK x Af(K), the balance equation ([6]) leads to the following integral equation 
for the surface elevation spectrum in steady situations, 

S nl [M(k)] + (eK x - ±a x e- a2 / 2 * 2 A M(k) = 0, (8) 



where the constant b — ^-p — is the duration of the breakin g pr ocess in units of r. There are two unknown param- 
eters, a and b. Contrary to similar equations in the literature [151 ] . the breaking- wave loss term is not proportional to 
a power of the amplitude of the fluctuations, it depends transcendentally on this amplitude, because wave-breaking 
is not found at any algebraic order in the amplitude expansion of the kinetic equations. Although g and U have been 
scaled out, the small dimensionless parameter e remains, in particular because of the term depending transcendentally 
on the unknown amplitude. This integral equation yields the scaling laws for the amplitude of the fluctuations by 
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noticing that the first term has "conservation laws". This is a familiar property of Boltzmann-type equations, which 
writes J dKF(K)S n i [A/"(K)] = with F(K) = 1, K or Q(K). Therefore one finds three relations to be satisfied by 
the steady spectrum : 



This yields the same scaling relation as derived before, if one assumes that the integration over K and the multiplication 
by F(.) change the scaling in the same way in the two terms of equation ©, the one representing the input by the 
instability (proportional to ek x U) and the one representing the loss by wave-breaking, proportional to the exponential. 
Notice too that, thanks to the prefactor a x in the dissipation part, the dissipation becomes more efficient as the 
amplitude of the spectrum grows. This makes likely impossible a steady (unphysical) spectrum spreading without 
decay to infinitely large wave numbers. 

To summarize, we have introduced the idea of a physical small parameter in the wind-sea interaction, the ratio of the 
mass densities of the two fluids. Thanks to that, we have been able to derive scaling laws for the physically observable 
quantities like the mean square slope of the sea surface. This gives a basis for the derivation of an Hassclmann-type 
equation for the steady state spectrum, that would be valid not only for the case of wealkly nonlinear turbulence, where 
it is generally adressed, but even for storms with very large wind velocities. We propose a fully explicit mathematical 
model, equation JS}, for the steady spectrum of surface elevation perturbed by a constant wind. This model is valid 
for large wind speeds, where dissipation is mostly due to wave-breaking, although at low speeds the wave-breaking 
is stopped (before the overhanging of the sea surface) by capillarity effects. Indeed all this remains to be confronted 
with experimental results alhough the solution of the spectral equation remains to be studied in details in the limit e 
small. 
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